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I.  Introduction 


In  a  recent  survey,  Parker  and  Kasemir  (1980)  address  the  need  for  reliable  and  inexpensive  ! 

airborne  I  ightning  warning  and  avoidance  systems.  They  survey  a  rnrtoer  of  techniques  currently 
in  use  for  lightning  research  and  their  potential  suitability  for  aircraft  systems.  Anong  those 
techniques  considered  for  warning  of  distant  storms  were  a  nuttier  of  sferics  detectors.  "Sferic" 

i 

has  been  used  to  cover  all  aspects  of  lightning  RF  radiation,  though  its  usage  usually  is  j 

limited  to  the  lower  frequency  (<10  hHz)  radiation.  Wtile  this  radiation  is  more  energetic  < 

than  the  higher  frequency  (l-F  and  WF)  radiation,  recent  advances  in  our  understanding  of  the 
VTF  radiation  from  lightning  suggest  that  it  might  be  feasible  to  use  the  VFF  radiation  in  the  i 

i 

development  of  an  airborne  warning  system 

The  present  study  was  undertaken  to  assess  the  feasibility  of  developing  an  airborne  light¬ 
ning  warning  system  which  detects  the  location  of  VFF  radiation  from  lightning.  Initially  a  I 

prototype  device  was  also  to  have  been  constructed,  so  the  study  focused  on  particular  target 
aircraft  which  were  to  be  available  for  testing  the  prototype.  Wien  funding  for  the  prototype 
was  cancel  led,  the  study  shifted  to  an  examination  of  a  generic  system  for  any  aircraft. 

The  phenomenological  characteristics  of  VFF  radiation  from  I  ightning  provide  definite  advan¬ 
tages  for  an  airborne  system  These  characteristics  are  examined  in  section  II.  First, 
relevant  results  f ran  studies  of  the  more  familiar  optical,  acoustic  and  electrical  emanations 
fron  lightning  are  discussed.  The  unique  characteristics  of  the  WF  radiation  are  then  presen¬ 
ted  with  arguments  for  their  use  with  an  airborne  syston  Section  II  concludes  with  a 
presentation  of  WF  lightning  parameters  Wtich  vrould  be  relevant  to  the  design  of  an  airborne 
system. 

The  other  factors  Wtich  trust  be  considered  for  the  feasibility  are  presented  in  two  catego¬ 
ries.  Section  III  considers  factors  associated  wi th  the  propagation  of  the  WF  lightning  radia¬ 
tion  from  its  source  to  the  aircraft.  Such  factors  as  signal  strengths,  dynamic  range,  interfer¬ 
ence  and  reflections  are  considered  Section  IV  considers  the  technological  aspects  of  a  system 
at  the  aircraft  itself.  The  suitability  for  aircraft  placement  of  existing  ground  techniques  is 


examined  for  both  azinuth  and  range  determinations.  Each  of  these  determinations  presents  tech¬ 
nical  prob  I  ems  which  mist  be  addressed  for  a  syston  It  is  shown  that  nuch  flexibility  exists  in 
the  choice  of  parameters  for  azinuth  determination.  The  available  possibilities  for  range  deter¬ 
mination  are  not  as  premising 

The  target  performance  for  this  study  of  an  airborne  system  was  left  open  to  our  determina¬ 
tion  We  set  a  desired  range  of  coverage  at  --150  Ian  For  a  craft  near  the  speed  of  sound  this 
translates  to  a  time  frame  of  roughly  10  minutes.  Angular  coverage  is  discussed  in  the  section 
on  azinuth  determination.  Nothing  was  provided  to  indicate  whether  this  system  would  be  used  in 
friendly  territory  only  with  possible  ground  support  or  whether  it  would  be  intended  for  use 
during  cortoat.  It  was  also  not  stated  whether  or  not  the  system  mist  be  purely  passive.  In  the 
conclusion  same  of  these  issues  wi  I  I  be  addressed  in  I  ight  of  vrtiat  we  bel  ieve  can  be  done  with  an 
airborne  RF  lightning  detector. 


II.  Phenomenological  Considerations 


A.  General  Considerations 

Lightning  has  been  defined  as  a  transient,  high-current  electric  discharge  vrtiose  path  length 
is  general  ly  measured  in  kilometers  (linan,  1969).  The  most  cornnon  occurrence  of  lightning,  and 
the  only  of  interest  for  this  report,  is  that  associated  with  thunderstorms.  Such  lightning  is 
usually  distinguished  between  that  which  takes  place  entirely  within  the  cloud  (intracloud  or 
cloud  discharges)  and  that  which  canes  in  contact  with  the  ground  (cloud-to-ground  or  ground 
discharges).  This  distinction  is  somewhat  misleading  in  that  all  discharges  have  significant 
activity  within  the  cloud.  Furthermore,  many  discharges,  especial  ly  those  frammature  and  carrplex 
storms,  appear  as  hybrids  or  cortiinations  of  cloud  and  ground  discharges. 

Although  the  most  frequently  occurring  form  of  lightning  is  that  within  the  cloud,  the 
greater  part  of  the  lightning  literature  concerns  the  ground  discharge,  especially  return 
strokes.  (An  entire  discharge  is  usually  referred  to  as  a  flash,  lasting  about  a  half-second; 
the  highly  liminous  components  of  discharges  fran  ground  to  cloud  are  referred  to  as  strokes, 
or  return  strokes,  lasting  on  the  order  of  milliseconds.)  These  distinctions  are  important 
because  the  incloud  processes  are  obviously  of  most  importance  frcm  the  standpoint  of  an  air¬ 
borne  platforms  Lhfortunate ly  much  of  the  work  on  airborne  systems  has  also  concentrated  on 
the  return  strokes  (e.&  Rustan  et  al,  1983).  Although  they  are  connected  in  the  overall 
structure  of  the  flash,  the  intent  of  this  discussion  on  I ightning  phenomenology  is  to  emphasize 
what  is  known  about  what  is  happening  in  the  cloud  where  an  aircraft  would  be  located  as 
opposed  to  the  point  vtftere  a  discharge  contacts  the  ground. 

Significant  advances  have  been  made  on  the  understanding  of  incloud  processes  within  the 
last  few  years.  Uhfortunate ly  much  work  remains  to  be  done  before  a  carplete  picture  wi  I  I 
anerge. 

Cptical  and  photographic  techniques  which  were  very  successful  in  the  determination  of 
processes  occurring  in  ground  discharges  have  been  of  little  use  in  studying  cloud  discharges, 
for  obvious  reasons.  During  a  cloud  discharge  a  continuous  low  luminosity  is  seen  with  several 
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relatively  bright  pulses  of  less  than  1  ms  duration  super  inposed.  Any  spatial  structure  can  not 
be  determined  due  to  scattering. 

Other  techniques  which  have  been  successful  in  determining  incloud  processes  are:  acoustic 
reconstruction  of  lightning  channels;  analysis  of  the  radar  returns  fran  lightning  channels 
within  the  cloud;  modelling  of  charge  transfers  from  the  analysis  of  charges  in  the  electric 
field  at  nunerous  ground  stations;  and  analysis  of  the  VFF  radiation  fran  lightning  to  deter¬ 
mine  the  location  of  sources  and  their  movement.  With  the  exception  of  some  earlier  electric 
field  studies,  real  progress  using  these  techniques  has  occurred  only  within  the  last  ten 
years.  All  studies  using  these  techniques  have  involved  a  small  nurber  of  flashes.  Usually 
little  documentation  has  been  available  on  the  thunderstorm  environment  producing  the  flashes. 
It  is  still  somewhat  premature  to  draw  general  conclusions  about  incloud  I  ightn  I  ng  processes 
and  their  connection  to  and  range  of  variability  with  respect  to  the  thunderstorm  environment. 
Nevertheless,  the  following  results  from  these  techniques  provide  the  phenomenological  back¬ 
ground  from  which  the  feasibility  of  an  airborne  platform  is  discussed. 

1.  Acoustic  Reconstruction.  Detailed  analyses  using  the  acoustic  technique  of  Few  (1970) 
have  been  presented  by  Teer  (1970)  and  MacGorman  (1978).  In  this  technique,  the  direction  of 
arrival  of  a  thunder  inpulse  at  an  array  of  microphones  is  determined  by  the  propagation  time 
between  microphones,  measured  to  an  accuracy  of  1  ns  using  cross-correlation  analysis.  The 
acoustic  ray  defined  by  this  direction  is  then  retraced  through  the  atmosphere  to  its  source. 
Range  to  the  source  is  determined  fran  the  propagation  time  between  occurrence  of  the  light¬ 
ning,  determined  by  an  electric  field  change  meter,  and  arrival  of  the  thunder  inpulse  at 
the  array.  Source  locations  are  corputed  from  all  processes  in  a  lightning  flash  that 
produce  thunder  at  the  array.  However,  since  acoustic  propagation  time  is  long  compared  to 
the  duration  of  lightning,  the  acoustic  analysis  cannot  separate  the  structure  of  individual 
processes  in  a  flash,  but  integrates  the  structure  from  a  I  1  processes. 

The  reconstructed  lightning  structure  was  showm  to  be  predominant  ly  horizontal,  usually 
nruch  larger  than  any  vertical  extent.  In  addition  the  flashes  for  the  10-25  minute  intervals  of 


study  tended  to  occur  in  one  or  two  layers  at  altitudes  above  the  0°  isotherm  In  same 
storm  the  horizontal  extent  of  the  flashes  increased  as  the  storm  progressed  fron  a  typical 
10-15  km  to  greater  than  20  km  The  increase  was  attributed  to  the  growth  of  new,  electrical¬ 
ly  active  cells  near  the  previous  activity.  It  should  be  noted  that  these  studies  were  conduc¬ 
ted  in  large,  usual  ly  mature  thunderstorms.  The  question  of  whether  all  lightning  from 
all  phases  of  a  thunderstorm  is  predominantly  horizontal  remains  open. 

2.  Radar  Observations.  Radar  studies  of  lightning  have  not  been  concerned  so  mich  with 
locating  lightning  within  the  storm  structure  as  with  determining  the  actual  properties  of  the 
lightning  channel.  Mazur  and  Rust  (1983)  present  some  data  concerning  lightning  propagation 
using  V1-F  and  L  band  radars.  They  separate  I  ightning  flashes  frcm  squal  I  I  ine  storms  in  Cklahona 
into  two  categories:  one  with  echoes  having  projected  lengths  onto  a  horizontal  plane  of  <20  km 
and  another  with  lengths  of  >20  km  These  lengths  are  determined  from  the  maximm  extent  of 
radar  return  from  the  ionized  channel  along  the  radar  beam  direction.  Their  study  shows  that 
the  largest  lightning  density  tends  to  be  near  the  leading  edge  of  the  precipitation  cores  in 
developing  cells.  As  a  cell  in  the  squall  I  ine  develops,  the  total  I  ightning  dens  i  ty  increases. 
Long  discharges  develop,  but  the  shorter  ones  predominate.  In  contrast,  as  the  cel  I  dissipates, 
short  flashes  diminish  or  cease  and  the  long  flashes  dominate  the  activity.  It  must  be 
errphasized  that  these  resul  ts  were  mainly  from  a  single  Oklahoma  squal  I  line  system 

3.  Electric  Field  Change  Studies.  Modern  observations  of  the  lectric  field  change  pro¬ 
duced  lightning  (Krehbiel,  et  al.,  1979;  Krehbiel,  1981)  have  contributed  greatly  to  our  under¬ 
standing  of  incloud  lightning  processes.  Measurements  are  made  of  the  electric  field  change  at  a 
nirrfeer  of  ground  stations  separated  by  a  few  kilometers  each.  These  are  then  fit  to  a  particular 
model  of  charge  change  vrtiich  would  give  the  measured  electric  field  changes.  Models  have  ranged 
frcm  a  sirrple  point  charge  model  involving  four  parameters  -  location  x,  y,  z,  and  magni  tude,  Q 
-  to  recent  work  (Liu  and  Krehbiel,  1984)  involving  an  advancing  streamer  model  with  seven 
unknown  parameters.  By  having  more  measurements  than  unknown  parameters,  the  goodness  of  the  fit 
between  data  and  model  can  be  examined  using  chi-squared  methods.  The  results  of  the  streamer 


model  indicate  that  the  initial  electric  field  change  of  the  intracloud  flashes  was  produced  by 
upward  motion  of  negative  charge  in  the  cloud  The  apparent  height  of  initiation  was  between  7.3 
and  8.3  ton  MSL  for  three  flashes  in  a  smal  I  storm  and  10  bn  for  an  energetic  flash  in  a  large 
storm  These  altitudes  correspond  to  enviroonental  torcperatures  of  -17  to  -24°C  and  -  38°C, 
respectively.  After  initiation  the  streamers  progressed  upward  with  a  speed  on  the  order  of  10^ 
m/s  or  greater.  In  all  cases  the  charge  transfer  at  the  beginning  of  the  flash  was  vertical  or 
nearly  vertical.  These  results  pertain  to  the  intial  15-30  ms  of  the  discharge;  after  that  time 
the  development  could  not  be  determined  with  the  models  used  Carparisons  with  WF  radiation 
analysis  (Proctor,  1976,  1981)  has  led  to  the  conclusion  that  lightning  begin  near  the  negative 
charge  region  of  a  thunderstorm  with  negative  streamer  propagation. 

Most  of  the  discharges  from  studies  discussed  previously  and  from  the  Proctor  studies  were 
horizontally  oriented  Liu  and  Krehbiel  attribute  discrepancies  to  the  fact  that  other 
studies  were  from  nature  or  dissipating  storms,  vtfiereas  the  flashes  fran  their  study 
occurred  in  graving  cells.  Krehbiel  (1981)  has  also  shown  several  extensive  horizontal 
discharges  from  a  dissipating  storm  system 

Krehbiel  et  al  (1984b)  present  a  study  of  the  initial  fifteen  flashes  in  a  small  Florida 
storm  The  electric  field  change  models  are  cmpared  with  radar  observations  of  precipitation 
development  and  velocities  in  the  storm  and  with  locations  of  VFF  sources  determined  by  the  KSC- 
LOVR  system  Figure  1  shows  a  plot  of  the  average  height  of  the  negative  and  positive  charge 
center  Jipole  model)  and  VfF  radiation  sources  for  the  fifteen  flashes.  The  ipper,  positive 
charge  centers  of  the  intracloud  lightning  are  located  above  15-20  dBZ  reflectivity  in  the  top 
of  the  cloud  and  increased  in  altitude  as  the  storm  grew  vertical  ly,  at  the  same  rate  (~8m/s) 
in  both  cases.  The  negative  charge  centers  remained  at  approximately  constant  elevation,  at 
tOTperatures  between  -10  and  -15°C  Two  cloud  to  ground  (GG)  flashes  removed  negative  charge 
from  essential  ly  the  same  altitude  and  volirne  of  the  storm  as  intracloud  (1C)  flashes  before  and 
after  then\  though  the  VFf  sources  indicate  that  the  CG  discharges  extended  downward  and  the 
1C  flashes  extended  upward  in  the  cloud 
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B.  VFF  Lightning  Phenomenology. 

In  order  to  utilize  the  VFF  radiation  from  lighting  in  a  location  system,  one  rrust  know 
both  the  terrporal  and  spatial  characteristics  of  that  radiation.  The  tenporal  characteristics 
are  important  because  of  differing  phenomenology  with  time  during  a  flash  and  because  the  terrpo¬ 
ral  characteristics  set  the  design  parameters  for  the  system  The  spatial  characteristics  are 
irrportant  for  two  reasons.  First,  the  sources  must  not  disrupt  the  measurement  of  their  loca¬ 
tion.  This  could  happen  if  there  were  numerous  widely  spaced  sources  radiating  during  a  single 
observation  period  (ie  "simultaneous"  sources)  or  if  the  source  positions  moved  significantly 
more  than  an  angular  resolution  element  during  the  observation  interval.  Second,  the  VFF 
sources  must  be  corre I ated  wi th  other  lightning  and  thunderstorm  phenomena  in  order  to  understand 
and  interpret  the  significance  of  their  location  Published  studies  using  the  ground  based 
techniques  of  a  hyperbolic  system  (Proctor,  1981,1983)  and  an  interferometric  syston  (Warwick 
et  al,  1979;  Hayenga  and  Warwick,  1981;  Ftoyenga,  1984)  provide  observations  fron  which  we  can 
infer  characteristics  of  the  spatial  structure  of  the  VFF  radiating  sources,  and  the  relation  of 
the  sources  to  other  lightning  and  thunderstorm  phenomena. 

These  observations  should  not  be  considered  comprehensive.  Much  work  remains  to  be  done  in 
this  area. 

1.  Tenporal  Structure  of  VFF  Radiation.  In  a  study  associated  wi  th  this  feasibility  study, 
Rhodes  (1985)  analyzed  seme  properties  of  lightning  radiation.  We  present  his  analysis  of  one 
flash  to  show  hcrw  VFF  radiation  evol  ves  in  time  and  how  it  is  related  to  other  I  ightning  phenome¬ 
na.  Lightning  was  observed  with  microsecond  resolution  at  four  RF  frequencies;  34,  368,  1414, 
2275  Mhz.  Supplemental  measurements  were  made  of  the  electric  field  change  and  the  optical 
radiation.  A  video  canera  also  captured  pictures  of  the  flashes. 

The  flash  presented  here  is  representative  of  the  flashes  studied.  Fkwever,  the  nirrber  of 
flashes  studied  was  too  small  to  rroke  an  adequate  assessment  of  the  range  of  variability  of  the 
parameters  discussed. 

The  flash  occurred  at  16  hours,  19  minutes,  23  seconds  on  September  18,  1982,  in  an  isolated 


ror  i  on  remains  fixed  in  elevation  whereas  the  positive  chaiy.e  region  moves  upward  at, 


thunderstorm  south  of  Socorro.  Figure  2  shews  two  frames  from  the  video  recording  of  this 
flaslv  The  time  insertion  unit  was  incorrectly  set  ten  seconds  fast.  The  time  to  thunder 
indicates  that  the  ground  channels  were  about  2  to  3  bn  from  the  observation  van. 

Figure  3  presents  an  overview  of  the  flash.  The  received  RF  powers,  optical  signal,  and 
electric  field  measurements  are  plotted  as  a  function  of  time  with  mil  lisecond  resolution.  The 
various  parts  of  the  flash  are  labelled  The  flash  begins  with  a  preliminary  breakdown  in  the 
cloud  This  in  most  easily  seen  in  the  IF  records  although  the  1414  and  2275  Miz  records  don't 
showmich  radiation;  probably  due  to  antenna  directionality.  Closer  examination  of  the  radiation 
during  the  preliminary  breakdown  phase  shows  that  it  consists  of  pulses  or  groups  of  pulses  of  a 
few  microseconds  separated  by  10  to  100  microseconds.  The  bursts  become  more  dense  with  time. 

The  preliminary  breakdown  leads  into  the  stepped  leader  as  the  incloud  streamer  turns  toward 
ground  The  apparent  increase  in  IF  power  is  due  to  the  integration  of  many  short  pulses,  now 
more  nunerous  than  in  the  preliminary  breakdown.  The  stepped  leader  stops  abruptly  at  the  first 
return  stroke,  indicated  by  the  large  anplitude  electric  field  change,  the  cessation  of  RF 
radiation  and  the  large  optical  signal.  This  first  return  stroke  Is  shown  in  Figure  2a. 

The  interstroke  period  shows  mostly  isolated  bursts  of  IF  radiation,  lasting  a  few  hundred 
microseconds  each.  Same  of  these  have  associated  optical  pulses  and  electric  field  signals. 
Such  electric  field  signals  have  been  named  K  changes. 

This  flash  is  somewhat  different  in  that  the  second  return  stroke  occurs  in  wAiat  was  a 
branch  of  the  first  return  stroke  channel,  as  shown  in  Figure  2b.  About  2  mi  I  I i seconds  after 
this  return  stroke  a  return  stroke  occurs  in  the  original  channel,  probably  initiated  by  the 
stroke  in  the  branch. 

During  the  times  around  the  return  stroke  the  RF  radiation  consists  of  continuous  radiation 
lasting  tens  to  hundreds  of  microseconds.  This  radiation  is  not  due  to  the  return  strokes  since 
it  usually  initiates  before  the  stroke.  The  radiation  may  continue  through  the  time  of  the 
stroke.  Figure  4  shows  data  near  the  time  of  the  second  stroke  with  higher  time  resolution. 
Since  this  was  a  new  channel  to  ground,  stepped  leader  radiation  can  also  be  seen  In  the  first 
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FIGURE  2.  VIDEO  PICTURES  OF  THE  SECOND  AND  THIRD  RETURN  STROKES  OF  FLASH  161923.  THE  TIME 
INSERTION  UNIT  WAS  INCORRECTLY  SET  10  SECONDS  FAST. 
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part  of  this  record. 


The  video  recording  shoved  that  the  third  return  stroke  also  occurred  in  the  branch  channel. 
The  expanded  time  plot  for  this  return  stroke  in  Figure  5  shows  a  corbination  of  short  inpulsive 
radiation  (-  a  few  ps)  and  quasi-continuous  radiation  (~  SO  to  a  few  hundred  us). 

The  last  thi  rty  mi  I  I  i  seconds  of  radiation  fron  this  flash  (cf.  figure  2)  are  fron  an  incloud 
streamer.  The  characteristics  of  this  radiation  are  similar  to  stepped  leader  radiation. 

These  Socorro  studies  basically  agree  with  French  studies  by  Le  Boulch  et  al  (1984)  and 
South  African  studies  by  Proctor  (1981). 

From  these  studies  it  has  been  concluded  that  lightning  RF  radiation  is  of  two  fundamental 
types.  The  first  consists  of  inpulses  lasting  a  microsecond  or  less  which  came  in  groups  lasting 
on  the  order  of  20-50  us.  Such  groups  are  separated  by  intervals  of  the  same  order  of  magnitude 
in  which  last  from  a  few  mi  I  I  i  seconds  to  several  tens  of  milliseconds.  The  second  type  of 
radiation  is  of  higher  amplitude  and  consists  of  essentially  continuous  radiation  for  periods 
from  tens  of  microseconds  up  to  one  mi  I  I  i  second  and  occurs  in  isolation  during  the  flash. 

Figure  6  shows  measured  pulse  durations  for  the  flash  previously  presented  at  four  differ¬ 
ent  frequencies.  The  durations  were  determined  fron  the  time  interval  the  detected  signal  was 
above  a  specified  threshold,  indicated  on  the  figure.  The  values,  in  <fh\  are  received  power  at 
receiver  input,  neglecting  antenna  gain  and  source  polarization  and  distance.  Only  the  first 
fifteen  time  bins  contained  data,  making  a  statistical  analysis  rather  difficult.  A  scattering 
of  pulse  duration  from  50  to  400  microseconds  was  measured,  indicative  of  the  second  type  of 
radiation.  There  is  not  a  significant  variation  in  pulse  widths  as  a  function  either  of  frequency 
or  threshold  level  other  than  the  clear  indication  that  lightning  radiated  strength  declines  at 
shorter  wavelengths. 

2.  \4-F  Spectrin  of  Lightning.  Published  reviews  of  measurements  of  the  spectrin  of  VFF 
radiation  from  lightning  are  similar  to  that  sham  in  Figure  7  from  Oh  (1969).  The  spectrum 
falls  off  at  1/f  or  greater  with  increasing  frequency.  Usually  the  only  attempt  to  accomodate 
the  tanporal  structure  of  the  radiation  is  by  normalization  of  the  various  receiver  bandwidths 


Figure  3.  Overview  of  flash  161923.  The  electric  field  change,  optical 
signal  and  detected  outputs  from  four  receivers  are  plotted.  The  2275 
MHz  signal  is  inverted.  The  preliminary  breakdown  begins  the  flash, 
leading  into  the  stepped  leader,  shown  by  the  dip  in  the  electric  field 
and  strong  RF  radiation.  The  large  discontinuities  in  the  electric  field 
are  due  to  return  strokes.  (From  Rhodes  1985.) 
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Figure  U.  Higher  resolution  plot  of  the  time  around  the  second 
return  stroke  of  flash  161923.  Stepped  leader  radiation  leading 
to  the  return  stroke  can  he  seen  in  the  3**  and  368  MHz  channels. 
The  optical  signal  saturated  at  the  return  stroke.  (From  Rhodes, 
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Figure  5-  Higher  resolution  plot  of  the  time  around  the  third  return  stroke  of  161923. 
Both  the  short,  impulsive  radiation  and  the  longer  continuous  radiation  can  be  seen 
(From  Rhodes,  1985. ) 


Figure  6.  Histograms  of  the  duration  of  pulses  for  radiation  at  four  different  frequencies 
Durations  were  determined  from  the  time  the  detected  output  stayed  above  the  given  signal 
level  threshold.  The  1  microsecond  bin  contains  all  pulses  less  than  one  microsecond  in 
duration. 


Figure  7.  Plot  of  the  spectrum  of  lightning  radiation  derived  by 
0h(l969)  from  various  observations.  As  pointed  out  in  the  text  this 
plot  is  somewhat  misleading. 
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«4)ich  are  sufficiently  narrow  to  integrate  mich  of  the  signal.  In  the  spectrin  shown  in  Figure  7 
the  points  from  Brook  and  Kitagawa  (1964)  and  Hewitt  (1957)  are  the  dominant  points  defining  the 
slope  of  the  spectrun  However,  the  levels  used  are  not  received  powers  but  actual  ly  the  minitrun 
detectable  signal  characteristics  of  their  receivers.  Because  the  spectrun  of  lightning  at  VhF 
is  so  poorly  specified  in  the  existing  literature,  we  attenpted  a  new  set  of  measurements  in  the 
experiment  by  Rhodes  (1985)  quoted  previously.  Although  attmpts  were  made  at  estimating  the 
received  power  spectrun  of  the  various  phases  of  the  lightning  flash,  the  results  were  of  limited 
value.  Large  uncertainties  in  the  receiver  calibrations  and  antenna  gains  together  with  a  lack 
of  knowledge  concerning  source  locations  and  polarization  made  accurate  estimates  of  power 
levels  inpossible. 

It  is  clear  that  there  are  still  no  re  1 1  able  measurements  of  the  power  spectrun  of  lightning 
at  VFF  frequencies.  These  measurements  are  somewhat  difficult  and  expensive  to  perform  How¬ 
ever,  srfien  design  parameters  are  discussed  in  Section  III,  we  wi  II  show  that  a  functional  system 
could  be  designed  without  these  measurements.  The  fine-tuning  of  such  parameters  could  be 
acconplished  during  the  testing  of  a  prototype  system 

3.  Locations  and  Movement  of  VHF  Sources.  It  is  inportant  to  have  some  understanding  of 
the  locations  and  movement  of  the  sources  of  VhF  radiation  from  lightning  for  two  reasons. 
First,  there  rmy  be  characteristics  of  their  locations  which  could  be  utilized  to  reduce  the 
complexity  or  increase  the  performance  of  the  location  system  If,  for  exanple,  al  I  the  sources 
from  the  flash,  or  a  recognizable  portion  of  it,  occurred  at  the  same  altitude,  the  need  for 
measuring  elevation  angles  could  be  reduced.  Second,  the  sources  should  not  move  sufficiently 
far  or  fast  as  to  confuse  the  measurement  of  their  positioa  In  order  to  obtain  an  accurate 
location  determination  we  will  need  to  average  over  a  nurrfeer  of  inpulses  spaced  in  time.  This 
incoherent  average  will  increase  the  location  accuracy  as  the  square  root  of  the  nurirer  of 
irrpulses  averaged,  assuning  all  inpulses  originate  in  nearly  the  same  place.  If  they  don't, 
the  location  accuracy  wi  I  I  not  go  up  as  the  square  root  and  nay  even  decline. 


\AF  radiating  lightning  sources  have  been  located  by  three  different  techniques,  all  based 
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on  time  differences  of  arrival  between  spatially  separated  antennas.  The  hyperbolic  system  uses 
antennas  separated  by  kilometers.  Equal  time  differences  of  arrival  are  situated  on  hyperbolic 
surfaces,  the  intersections  of  wrfiich  provide  the  source  location.  By  using  antennas  spaced 
meters  apart,  the  hyperbolic  surfaces  becane  cones  (i.e.,  direction  cosines),  sinplifying  the 
caiputational  geometry.  Each  of  these  time-di fference-of-arr I val  systems  requires  recognizable 
anplitude  structure  in  order  to  make  a  measursnent.  The  interferometric  method  eliminates  the 
need  for  recognizable  snplitude  structures  by  measuring  the  phase  difference  between  signals  at 
the  spaced  antennas. 

Two  hyperbolic  systems  have  been  used  for  Wf  lightning  studies,  the  Lightning  Detection  and 
Ranging  (LDAR)  system  at  KSC  and  a  system  bui  1 1  by  Proctor  in  South  Afr  ica.  The  LDAR  had 
provided  a  few  data  points  per  flash  in  real  time  al  lowing  for  an  overview  of  many  flashes  in  a 
storm  whereas  the  South  African  sysum  has  been  used  to  provide  high  detail  of  a  few  flashes. 
Detailed  analysis  of  the  LDAR  data  (festan,  1979)  has  been  discredited  (Krehbiel  et  al,  1984a). 

Figure  8a  shows  LDAR  data  fran  an  entire  storm  showing  height  as  a  function  of  time.  Indi¬ 
vidual  flashes  cannot  be  resolved.  The  LDAR  sources  were  produced  primarily  by  intracloud 
discharges  and  tended  to  be  biased  toward  the  upper  or  positive  end  of  the  discharge.  Most  of 
the  sources  were  above  8  km  (MSL)  altitude  and  below  12  Ion  (MSL)  altitude.  The  upward  sloping 
radiation  patterns  were  correlated  with  the  occurence  of  intracloud  lightning  in  individual, 
growing  convective  cel  Is.  This  can  be  seen  in  Figure  8b  vrftich  is  a  plan  view  of  the  radiation 
sources  in  during  a  ten  minute  interval.  Three  cells  are  evident. 

Very  few  points  shown  in  Figure  8a  are  from  channels  to  ground.  This  is  in  agreement  with 
the  tenporal  data  presented  in  the  last  section  in  which  the  channel  to  ground  radiated  only 
during  a  srml  I  portion  of  the  overal  I  VFF  radiation. 

A  consistent  feature  which  merges  frem  the  analysis  of  a  few  more  storms  is  the  8  Ion  (MSL) 
threshold  for  the  LDAR  points.  Cells  did  not  produce  lightning  until  significant  radar  echoes 
extended  to  8  tan  altitude  in  the  cloud  and  the  negative  charge  center  for  the  lightning  appear  to 
remain  fixed  at  a  constant  altitude  (about  7  or  -15  degrees  Q  as  the  cell  grew  (cf  Figure 


1  )•  Although  we  can  set  a  lower  threshold  for  the  predoninant  VFF  radiation,  a  spread  (in 
elevation)  of  four  to  six  ki lameters  sti I  I  exists. 

In  the  previous  section  two  types  of  VTf  radiation  were  noted.  Proctor  (1981),  using  a 
hyperbolic  systor\  and  Hayenga  and  Warwick  (1981),  using  an  inter feronetr ic  systcn\  have  deter¬ 
mined  propagation  speeds  for  the  two  types. 

Hie  radiation  from  the  short  inpulsive  radiation  extends  at  about  105  rtVs.  Figure  9 
shows  this  drifting  motion  of  nunerous  sources.  The  diameter  of  the  drifting  region  is  on  the 
order  of  600  meters,  tote  the  horizontal  motion  of  the  preliminary  breakdown)  from  0-30  ms  and 
the  downward  motion  of  the  stepped  leader  at  30-44  ms. 

The  longer  duration  radiation  has  been  labeled  Q  noise  by  Proctor  and  fast  bursts  by 
Hayenga.  During  the  few  hundred  microseconds  of  radiation  these  bursts  travel  at  speeds  of  about 
2.5x107  nVs,  two  orders  of  magnitude  faster  than  the  short  inrpulsive  sources.  Such  a  burst  is 
shown  in  Figure  1  Ol  Same  care  must  be  exercised  In  determining  the  speed  since  it  involves 
differentiating  noisy  data.  The  point-to-point  length  of  the  path  shown  is  17  bn,  yielding  a 
speed  of  6x107  nVs.  Along  the  smooth  drawn  in  path,  representing  the  general  motion,  the 
length  is-*)  bn  giving  a  speed  of  2x107nrVs,  vrtiich  we  consider  more  representative. 

C.  Surma  ry 

The  results  from  these  studies  which  are  relevant  to  the  feasibility  of  an  airborne  system 
can  be  suimarized  as  fol  lows: 

1.  Lightning  appears  to  begin  with  a  negative  streamer  near  the  negatively-charged  region 
of  the  storm  This  region  is  around  O0-!^  temperature.  Its  actual  altitude  appears  to  be 
storm-dependent.  During  this  phase  VFF  sources  radiation  canes  in  short  pulses  whose  sources 
move  at  10^  nV$. 

2.  Ground  discharges  appear  to  move  dcumward  from  this  region.  The  more  cannon  intracloud 
discharges  appear  to  move  upward  from  this  region. 

3.  During  many  phases  of  storms  lightning  rmy  have  extensive  horizontal  extent  within  the 
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Figure  8.  VHF  radiation  source  data  for  a  storm  at  K.'r  from  the 

real-time  .TAP  system,  a)  Time  height  nlot  shoving  unward  trend 

of  mints  with  nil  growth .  Note  that  most  points  lie  between  8 
and  12  km  MSL.  h)  Plan  view  of  ten  minutes  of  data  showing  the 

points  grouped  in  cells.  ( Prom  Krehbiel  et  al,  198ka.) 
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4*  VrF  radiation  also  has  bursts  of  hundreds  of  microseconds  duration  vtfiose  sources  move  at 
speeds  in  excess  of  10^  nVs- 

5.  Most  VFF  sources  are  at  elevations  >8  km  MSL  with  a  spread  of  sane  4  ki  laTBters.  It  nust 
be  emphasized  that  these  results  are  still  new  and  not  yet  comprehensive.  It  is  not  at  all  clear 
what  relation  lightning  activity  at  a  given  point  in  time  may  have  with  that  10  minutes  later, 
for  example.  It  can  be  concluded  that  the  negatively-charged  region  may  remain  intact.  However, 
the  entire  region  enccrrpassed  by  lightning  may  change  dramatically  over  such  an  interval,  in¬ 
creasing  (or  even  decreasing)  by  kilaneters  in  both  horizontal  and  vertical  extent.  If  it  were 
necessary  to  fly  a  craft  within  2-5  Ion  of  an  active  lightning  region  it  would  be  necessary  to 
monitor  the  activity  on  approach  rather  than  rely  on  data  which  would  beccme  quickly  outdated. 

Wiile  the  results  presented  give  same  indication  of  where  lightning  is  occurring  within  a 
storn\  there  is  no  data  to  indicate  where  it  will  not  occur.  Neither  is  there  enough  data  fron 
these  studies  to  develop  statistics  concerning  the  probability  of  a  strike  to  a  craft  as  function 
of  tenperature  or  possibly  elevation.  Studies  by  N4SA  using  an  F106B  (eg.  Fisher  and  Pinner, 
1983)  may  be  more  helpful  in  this  regard. 

A  very  important  aspect  which  onerges  from  the  phenaneno I og i ca I  considerations  is  the  need 
to  detect  all  discharges,  not  just  the  return  strokes  of  ground  discharges.  Many  storms  show 
long  interval  s  with  very  active  incloud  flashes  and  very  few  ground  strokes.  No  good  correlation 
exists  between  the  location  of  ground  strokes  and  incloud  activity.  Aside  fron  the  possible 
errors  in  actual  operation  of  the  Stormscope  instriment  presented  by  Parker  and  Kasemir  (1982), 
phenomenological  considerations  indicate  that  the  Stormscope  would  not  adequately  detect  the 
extent  of  lightning  act i v i ty  wi thin  the  cloud  and  might  have  a  lew  probabi I ity  of  detection  for 
sore  storms. 

The  Stormscope  uses  the  vertically-polarized  50  kHz  radiation  from  lightning.  Such  radia¬ 
tion  comes  primarily  fron  return  strokes.  4ny  cloud  radiation  would  not  follcw  the  model  used 
for  range  detection  which  is  based  on  return  storkes.  It  is  fortuitous  that  the  errors  in  the 
syston  lead  to  a  larger  than  actually  detected  region  of  avoidance  providing  a  conservative 


ts  showing  a  spread  of  about  600  meters  in  the  points,  b)  Millisecond  averages  of  the 
lotion  of  2x10'  m/s.  (From  Havenga  and  Warwick,  1981 . ) 


warning  of  severe  weather  as  indicated  by  Baun  and  Seymour  (1979).  Amore  reliable  system  should 
detect  the  inctoud  hazard  as  it  develops  In  the  path  of  the  craft.  The  survey  of  the  VhF 
phenomenology  of  lightning  shows  why  it  is  the  preferred  activity  to  detect. 
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III.  Propagation  Conti  derat  ions 

A.  General  Considerations 

The  RF  signal  from  a  I  ightning  source  can  be  treated  as  any  other  signal  in  terms  of 
propagation  to  the  aircraft,  with  a  fo*  notable  exceptions.  Most  of  the  fundamental  work  on 
propagation  at  the  frequencies  of  interest  ms  performed  many  years  ago  as  reported  by  Kerr 
(1951)  and  Reed  and  Russell  (1953).  It  is  not  the  intention  to  rehash  any  of  that  work  here,  but 
rather  to  present  aspects  of  it  which  are  relevant  to  this  study  and  refer  the  reader  to  those 
texts  for  more  detai  I  s. 

A  single  propagation  path  consisting  of  a  direct  path  and  a  ground  reflection  path  between 
the  lightning  source  and  the  aircraft  is  assured  for  this  study.  This  path  is  shown  in  Figure 
11.  he  distance  between  the  source  and  the  craft  ranges  from  10  to  150  tar\  the  aircraft  altitude 
varies  from  5000  ft  to  30,000  ft  above  ground  and  the  lightning  source  varies  from  2  to  10  km 
above  ground.  A  plane  geometry  (flat  earth)  is  also  suitable  for  this  study  although  corrections 
might  be  needed  in  a  prototype  phase  to  account  for  spherical  geometry. 

Appendix  I  shows  the  elevation  angle  of  arrival  for  sources  within  the  selected  geometry. 
Plots  were  rmde  for  each  aircraft  altitude  fron  5000  ft  to  30,000  ft  in  5000  ft  increments.  On 
each  plot  the  elevation  angle  of  arrival  for  the  direct  and  ground  ref  I ected  waves  is  plotted  as 
a  function  of  range  from  10  to  150  km  for  I  ightning  sources  at  2,  4,  6,  8,  and  10  km  above 
ground.  The  range  of  elevation  angles  plotted  is  +_  20°.  These  plots  demonstrate  why  the 
ref  I  ected  wave  nust  be  considered  At  higher  aircraft  altitudes  direct  waves  are  received  at  the 
same  angles  at  which  reflected  waves  are  received  when  the  plane  is  at  lower  altitude.  The 
possibly  of  having  a  directional  antenna  which  is  steered  as  a  function  of  aircraft  altitude  was 
not  considered  as  it  would  unduly  carpi  icate  the  system 

Anomalous  propagation  such  as  ducting,  troposcatter,  diffraction,  or  rruitiple  reflections  in 
not  considered  Any  such  propagation  would  cause  an  Incorrect  position  determination  and  would 
degrade  the  performance  of  the  system  Most  anonalous  propagation  would  affect  the  apparent 
elevation  of  sources  with  the  exception  of  nul  tiple  reflections  (such  as  fron  mouta  ins  ides)  which 


pure  11.  System  peometry  assumed  for  this  study.  It  is  assumed  two  waves  will  be  received  at 
ie  aircraft;  a  direct  wave,  and  a  specularly  reflected  pround  wave.  Parameters  for  the  system 
e  discussed  in  the  text. 


would  affect  azlnuth  determination.  This  propagation  is  difficult  to  characterize  with  respect 
to  the  general  operation  of  this  syston  and  should  be  left  for  the  prototype  stage  in  order  to 
assess  the  inpact  on  performance  with  actual  data. 


Using  the  assured  propagation  path  together  with  the  phenomenological  considerations  of  the 
preceeding  sections  it  is  possible  to  parameterize  the  operating  frequency,  required  mininun 
detection  level  and  dynamic  range  required  of  the  syston  The  effects  of  the  ground  reflection 
path  can  also  be  assessed  Most  of  this  material  can  be  found  in  standard  references  (Kerr, 
1951;  Reed  and  ftissel,  1953;  and  Erst,  1984)  with  only  a  few  significant  differences. 

The  most  inportant  difference  to  consider  is  that  the  reference  texts  general  ly  deal  with  CW 
signals,  i.e.  signals  of  duration  long  compared  to  the  time  for  propagation;  whereas  lightning 
radiation  is  very  inpulsive.  Instead  of  interference  lobe  structures,  the  designer  of  this 
system  mist  consider  miltiple  pulses  from  various  regions  of  space  with  duration  and  duty  cycle 
such  that  they  may  or  may  not  be  received  sinul  taneous  ly  at  the  aircraft. 

Appendix  II  shows  plots  of  the  difference  in  path  length  between  direct  and  reflected  path 
for  the  geometries  of  interest.  Path  difference  and  equivalent  time  of  arrival  difference  are 
plotted  as  a  function  of  range  for  the  same  parameters  of  aircraft  elevation,  range,  and  source 
height  used  in  the  plots  in  Appendix  I.  A  caparison  of  the  time  differences  shown  and  the  pulse 
durations  and  intervals  between  pulses  given  in  the  previous  section  shows  that  the  systemwill 
receive  a  cent) i nit  ion  of  isolated  pulses  from  either  the  direct  of  reflected  waves  or  simulta¬ 
neous  radiation  from  both  sources  with  the  direct  wave  radiation  starting  before  the  similtaneous 
radiation  and  ref  lected  wave  continuing  on  afterward  Any  location  system  mist  be  able  to  handle 
any  of  these  situations  in  order  to  perform  reliably. 

B.  System  Design  Parameters 

The  three  design  parameters  to  consider  in  this  section  are  operating  frequency,  receiver 
sensitivity  or  mininun  detectable  signal  (M>5),  and  dynamic  range  of  the  system  Since  these 
parameters  wi  I  I  be  interdependent  in  this  systemwewill  consider  them  together. 

The  expected  signal  strength  at  the  receiving  antenna  will  allow  a  determination  of  the 
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necessary  sensitivity  of  the  system  Oily  the  direct  path  wi  I  I  be  considered  initially.  The 
power  received,  RP,  is  given  by 

RP  =  ERP  -  PL, 

where  ERP  is  the  effective  radiated  power  of  the  source  and  PL  is  the  path  loss,  a  I  I  measured  in 
dB.  Assuring  the  reference  standard  of  an  isotropic  receiving  antenna,  the  path  loss  is  given  as 
(Erst,  1984) 

PL  =  92.45  +  20  log  f  +  20  log  d, 

where  d  is  the  path  length  in  kiloneters  and  f  is  the  frequency  in  Q4z. 

Over  the  path  lengths  of  interest,  fran  10  km  to  150  tor\  the  path  loss  varies  by  23.5  <fi. 
The  data  fran  Rhodes  (1985)  indicate  a  variation  in  lightning  signals  of  less  than  40  dB.  We  can 
conclude  that  a  60  cfi  dynamic  range  should  be  adequate  for  the  lightning  system  A  60  cB  dynamic 
range  is  easily  achieved  with  standard  receiver  design  techniques. 

The  variation  in  path  loss  as  a  function  of  frequency  is  20  cfi  fran  30  kHz  to  300  kHz.  In 
the  absence  of  good  spectral  infomation  the  choice  of  operating  frequency  in  difficult  to 
specify,  but  a  change  of  about  20  dB  in  expected  ERP  fran  30  kHz  to  300  kHz  is  reasonable.  A 
system  at  30  kHz  would  need  to  be  about  40  dB  less  sensitive  than  one  at  300  MHz.  In  a  later 
section,  the  need  for  a  higher  frequency  of  operation  in  order  to  reduce  antenna  array  sizes  is 
stressed  as  a  trade-off  with  the  push  toward  lower  frequencies  to  reduce  requironents  on  system 
sensitivity.  In  a  prototype  system  design,  it  is  likely  that  neither  factor  will  define  the 
operating  frequency.  Any  syston  considered  in  this  report  will  require  a  bandwidth  of  one 
megahertz  or  greater,  thus  (raking  the  system  susceptab I e  to  interference,  from  both  transmitters 
on  the  ground  and  fran  other  systems  on  the  craft.  Freeing  the  syston  fran  interference  would  be 
the  strongest  factor  influencing  the  selection  of  operating  frequency.  It  should  be  noted  that 
Increasing  the  frequency  wi  1 1  require  a  lower  MX  in  order  to  receive  lightning  at  the  greater 
distances,  thus  again  increasing  susceptab) ti ty  to  interference. 


As  in  the  choice  of  operating  frequency  it  is  difficult  to  determine  the  required  mininun 
detectable  signal  (MS)  for  this  system  without  accurate  information  on  the  effective  radiated 
power  fron  lighning.  Wile  can  make  some  estirrate  from  the  w>rk  of  Rhodes  (1985).  His  system  at  34 
Mb  had  an  M3S  of  -85  (fin  while  the  syston  at  368  M-b  had  an  NOS  of  -100  (fin  Both  received 
lightning  signals  out  to  at  least  20  Ion  In  order  to  receive  lightning  out  to  a  direct  line  of 
sight  path  of  ISO  Ion  would  require  an  additonal  17.5  dB  of  sensitivity,  or  a  -  102.5  <fin  NOS  at 
30  Nttz  and  -117.5  (fin  MS  at  300  M-b.  Assiming  a  1  Mb  bandwidth,  the  30  M-b  system  could  have 
a  noise  figure  up  to  11.5  (fi  whereas  a  300  Mb  system  would  be  unattainable  without  very  compli¬ 
cated  electronics.  From  these  estimates,  it  can  be  concluded  that  the  operating  frequency  nay 
need  to  be  kept  below  100  Mb  to  avoid  syston  complexities  in  achieving  a  reasonable  MS  to 
receive  lightning  signals  out  to  150  bn  Better  definition  of  these  paranvters  would  have  to  be 
done  orpirical  ly  during  a  prototype  stage. 

During  the  prototype  phase  the  effect  of  ground  reflections  could  also  be  assessed.  While 
there  is  considerable  discussion  of  reflections  in  the  literature  (Reed  and  Russel  I,  1953;  Kerr, 
1951),  the  effects  from  the  ground  on  an  inpulsive,  rapidly  moving  source  are  not  clear.  While 
the  losses  on  ref  lection  may  be  sufficient  to  put  the  reflected  signal  below  detectable  levels, 


IV.  Technological  Considerations 

A.  General  Considerations 

In  the  previous  two  sections  we  have  examined  the  nature  of  the  WF  lightning  signal  as  it 
originates  at  the  source  and  the  effects  of  propagation  to  the  aircraft.  This  section  examines 
the  detection  of  that  signal  and  the  processing  of  the  signal  to  provide  information  on  its 
I  oca  t  ion.  The  determination  of  location  is  discussed  as  two  separate  determinations,  azimuth  and 
range.  The  azimuth  determination  covers  both  bearing  and  elevation  above  or  below  that  of  the 
craft.  Range  determination  covers  finding  the  distance  between  the  source  and  craft. 

The  discussion  of  location  determinations  is  made  from  a  systems  point  of  via*  Parameters 
such  as  accuracy,  angular  coverage,  etc.  are  discussed  without  reference  to  particular  hardware 
imp  lamentations.  Effects  such  as  errors  from  hardware  are  more  appropriate  for  a  prototype 
development  phase. 

B.  Azimuth  Determination 

1.  Expected  Resolution.  In  general,  direction  finding  (DF)  systems  fit  in  two  broad 
categories.  The  first  involves  comparisons  of  signal  amplitudes  received  with  kimm  antenna 
patterns  to  determine  the  direction  of  arrival  which  would  produce  those  amplitudes.  The 
crossed  loop  direction  finder  (QCF)  such  as  Stomscope  is  a  simple  example  of  such  a  system  To 
our  knowledge  no  antenna  pattern  systems  have  been  used  at  WF  for  lightning  location,  although 
there  are  current  examples  of  such  systors  at  WF  being  used  In  military  applications.  For 
satisfactory  operation  a  high  degree  of  antenna  pattern  accuracy  is  required,  along  with  predict¬ 
able  source  polarizations.  Carter  (1957)  has  shown  that  this  type  of  direction  finding  is  not 
feasible  for  aircraft  imp) emenation  at  hF,  especially  on  propeller  craft.  Most  reasons  given 
are  also  applicable  at  WF  and  for  those  reasons  and  the  lack  of  good  polarization  characteris¬ 
tics  of  WF  lightning  radiation,  such  systems  were  not  considered  feasible. 

The  second  broad  category  of  CF  systems  involves  the  correlation  of  signals  received  at 
spatially  separated  antennas.  Systems  correlating  the  amplitudes  are  usually  termed  Time- 
Difference-of-Arri val  (HA  or  TEA)  systems.  Those  using  a  phase  correlation  are  interferometers. 


Both  TO  and  interferometric  systems  have  been  used  for  lightning  location  at  \*F  frequencies  and 
have  been  proposed  for  aircraft  inp  lamentation  (Parker  and  Kasemir,  1982,  Section  V.  A  A  B). 
In  their  discussion,  Parker  and  Kasemir  present  the  expected  bearing  errors  for  the  two  systems, 
although  the  exanple  parameters  are  not  nearly  optlmrv 

For  a  time  of  arrival  system  the  angular  error  (A0)  is  related  to  the  time  resolution  for  a 
given  pulse  (At)  approximately  by 

o  60c  At 


where  c  is  the  speed  of  light  and  d  is  the  baseline  dimension  (antenna  separation).  Note  that 
this  is  independent  of  operating  frequency.  The  two  ways  to  decrease  angular  errors  are  by 
increasing  d,  which  is  difficult  to  do  in  an  aircraft  system,  or  to  decrease  At,  that  is 
increase  the  time  resolution.  Time  resolution  can  be  increased  by  an  increase  in  system  band¬ 
width.  This  will  lead  to  increased  interference  frcm  non-l  ightning  noise.  Moreover,  the  pulse 
rise  time  will  be  the  limiting  factor  in  deteimining  time  resolution.  The  lightning  signal  mist 
contain  short  rise-time  pulses  of  sufficient  energy  to  be  received  ig>  to  150  km  away  by  a  system 
with  sufficient  bandwidth  to  not  reduce  that  rise-time.  It  is  not  clear  that  such  pulses  occur 
within  lightning  radiation  with  any  regularity.  The  severity  of  retired  pulse  shape  and  the 
susceptibility  to  interference  are  sufficient  to  limit  the  feasibility  of  TCA  systems  for  azimuth 
determination. 

The  accuracy  of  angle  determination  in  an  interferometric  system  depends  in  four  parameters 
of  the  system;  operating  wavelength,  X  ;  baseline,  d;  system  bandwidth,  B;  and  the  integration  or 
averaging  time  interval,  T,  assuming  a  signal  is  present  during  the  entire  averaging 
period  The  expected  phase  measurement  error,  in  degrees,  is 

90  X  1 
77  3  72B-1  ’ 

as  derived  by  Hayenga  (1979).  Angular  errors  are  of  the  some  order,  but  are  dependent  on  actual 
direction  of  arrival  as  shown  by  Jacobs  and  Ralston  (1981). 


The  four  parameters  for  determining  the  accuracy  of  angular  measurement  for  an  interferome¬ 
tric  syston  are  not  only  interrelated  but  also  are  related  to  many  other  aspects  of  system 
design,  such  as  angular  coverage,  system  canp lex ity  and  system  size.  Mich  flexibility  exists  in 
the  definition  of  a  system  by  trading-off  various  aspects  of  the  system  Before  indicating 
choices  for  any  particular  system  we  will  spell  out  the  various  trade-offs  possible  with  these 
parameters. 

In  the  previous  sections,  we  presented  phenomenological  and  propagational  reasons  for  using 
a  lower  frequency  of  operation  for  a  location  system  For  an  inteferametr ic  system,  longer  X  re¬ 
quires  that  the  baseline,  d  also  be  increased  in  order  to  maintain  the  same  accuracy  and  requires 
that  antenna  size  increase  to  maintain  system  sensitivity.  Both  impact  the  size  of  the  system 
require  a  trade-off  with  aircraft  size. 

For  a  system  receiving  a  hemisphere  with  a  baseline  longer  than  X/2,  the  phase  measurement 
will  be  sibiguous,  since  phase  is  only  measured  modulo  2.  Jacobs  and  Ralston  (1961)  discuss 
methods  for  resolving  this  anbiguity  by  adding  antennas.  Arbiguity  resolution  roughly  double  the 
carpi  ex  ity  of  a  system  so  initial  ly  we  wi  1 1  consider  only  a  system  with  X/2  basel  ine  length.  At 
30  ME,  this  would  mean  a  baseline  of  5  meters,  reasonable  for  placement  on  most  aircraft  wings. 
Below  30  NHL,  aircraft  interactions  becane  inportant  (Carter,  1957). 

With  the  limitation  to  a  X/2  baseline,  the  error  equation  now  is  dependant  on  the  inverse 
square  root  of  the  product  of  B  and'T.  Precision  is  increased  by  maximizing  the  BT  product. 
Hayenga  (1979)  and  Richard  and  Affray  (1985)  have  shown  that  B  must  be  no  more  that  one  tenth 
the  operating  frequency  to  maintain  signal  coherency,  and  ideally  should  be  less.  In  the  VHF 
band,  B  would  range  frem  3  Mhz  at  30  Mhz  to  30  Mu  at  300  Mu.  However,  it  is  certain  that  for  a 
syston  to  operate  in  the  general  enviroonent  throughout  which  an  aircraft  would  fly  interfer¬ 
ence  fran  non-l  ightning  sources  will  limit  B  mich  further.  It  is  more  feasible  to  consider  1  Mhz 
as  a  rraximm  reliable  bandwidth.  Choice  of  bandwidth  and  operating  frequency  will  have  to  be 
chosen  careful  ly  during  prototype  development  to  avoid  interference. 

With  the  assumed  1  Mu  bandwidth  and  X/2  baseline  the  error  equation  is  dependent  only  on T. 


Representative  values  of  A6  for  various  values  of  Tare: 


T 

A8 

1  us 

21° 

10  ys 

6.5° 

100  us 

2.1° 

1  ms 

.65° 

At  150  ton  (worst  case)  range,  the  2.1°  resolution  obtained  with  a  100us  averaging  interval 
corresponds  to  a  region  5.5  ton  across.  Phenomenological  considerations  have  been  presented 
earlier  shoving  the  need  to  monitor  regions  of  activity  as  the  craft  approaches  to  note  any 
evolution  in  the  lightning  producing  regions.  By  the  time  the  craft  is  within  50  ton  of  the 
lightning  activity  the  2.1°  resolution  now  corresponds  to  a  region  less  than  2  ton  across.  This 
is  less  than  a  tolerable  distance  for  avoidance  of  lightning  within  clouds.  It  is  safe  to 
conclude  then  that  a  reliable  2°  resolution  in  aaimith  is  sufficient  for  the  operation  of  this 
system  and  would  be  attainable  with  a  100  us  averaging  time.  However,  the  error  equation  is  only 
valid  if  all  sources  occurr  ing  within  an  averaging  interval,  T,  arrive  from  the  same  angle^Afi. 
Miltiple  (stmil taneous)  sources  or  rapidly-moving  sources  will  invalidate  the  equation.  In 
research  inplementations  T,  has  been  kept  very  short  (-1-2  ys)  in  order  to  avoid  these  prob¬ 
lem.  Since  an  airborne  system  does  not  need  such  high  time  resolution,  we  can  use  other  methods 
to  deal  with  potential  problems  frcmnultiple  or  moving  sources. 

The  short,  inpulsive  radiation  whose  sources  move  atiMO^  nVs  will  not  present  a  major 
problem  The  speed  represents  only  a  movement  of  10  meters  in  100  microseconds,  with  the  spread 
of  sources  around  this  motion  an  order  of  magnitude  larger.  Ao  average  of  these  points  will 
locate  the  proper  region  of  space.  Since  the  duty  cycle  during  these  periods  of  radiation  is 
about  20%  (Hayenga  and  Warwick,  1981),  the  actual  interval  of  averaging  would  be  on  the  order  of 
500  microseconds.  The  occurence  of  the  pulses  could  be  determined  with  an  arrplitude  threshold 
criteria,  but  this  would  not  be  necessary  as  the  phase  measurements  of  the  noise  are  random  and 
would  cancel  in  the  average. 

The  fast-moving  («*10^  nVs)  sources  of  radiation  would  cover  1  ton  in  the  100  microsecond 


averaging  time.  A  single  data  point  at  the  average  position  would  represent  this  motion  and 
might  be  misleading  for  nearby  sources.  The  1  bn  source  extent  would  cover  more  than  the 
assured  2°  resolution  for  sources  closer  than  28  bn  aw ay. 

hfciltiple  sources  would  occur  for  two  reasons.  First,  there  is  the  possibility  of  radiation 
occurring  in  different  parts  of  the  sky  during  the  100  microsecond  averaging  time  (or  500  micro¬ 
seconds  if  the  duty  cycle  is  considered).  Second,  the  ground  wave  and  direct  wave  sources  wi  I  I 
carbine  to  form  three  distinct  sources;  the  direct  wave  alone,  the  ground  wave  alone,  and  the 
superposition  of  the  two  sources.  With  a  simple  interferometer,  it  is  not  possible  to  determine 
from  a  phase  measurement  the  miltiple  sources  giving  that  phase  measurement.  However,  it  is 
possible  to  determine  the  phase  measurement  vrtiich  will  result  fran  a  given  distribution  of 
sources.  The  details  are  presented  in  the  Appendix  to  Hayenga  and  Warwick  (1981).  It  is  shown 
that  for  two  sources  of  small  separation  (<10°),  the  result  of  these  sources  is  a  single  source 
location  between  the  two.  If  we  allow  smalt  deviations  in  source  position  but  reject  large 
deviations,  the  correct  region  of  lightning  occurence  can  be  identified. 

It  appears  feasible  to  overcame  potential  problems  fran  moving  or  nrultiple  sources  in  the 
fol  lowing  manner.  Tests  run  on  a  IMS  320  evaluation  board  and  experience  with  the  new  generation 
ground-based  interferometer  (Hayenga,  1983)  have  shawi  that  it  is  possible  to  process  interfero¬ 
meter  (Hayenga,  1983)  have  shown  that  it  is  possible  to  process  interferometer  data  in  a  pipe¬ 
lined  fashion  in  less  than  5  microseconds  per  location.  Assume  we  process  every  10  microseconds 
of  data.  The  expected  error  is  6.5°  per  data  point.  If  we  store  same  50  data  points  (ie  500 
microseconds  of  data)  we  can  perform  a  Randan  Strple  Consensus  (RAN5AC)  as  described  by  Fischler 
and  Bol  les  (1981)  to  determine  the  subset  of  those  data  points  which  came  from  the  same  region  of 
space.  The  algorithm  el  iminates  data  points  far  away  from  an  average  of  the  data.  Most  averag¬ 
ing  techniques  rely  on  the  assurption  (smoothing  assurption)  that  the  rnaxirrun  expected  deviation 
of  any  datun  fran  an  assumed  model  is  a  direct  function  of  the  size  of  the  data  set  and  that 
there  will  always  be  enough  good  data  to  smooth  out  gross  deviations.  As  this  nay  not  hold  for 
moving  or  nul  tiple  sources,  the  R*NSAC  algoritlvn  is  wel  I  suited  to  this  appl  i  cat  ion.  The  algo- 
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r i thru  has  three  unspecified  parameters;  the  error  tolerance  used  to  determine  whether  or  not  a 
point  should  be  accepted,  the  ntnt>er  of  subsets  which  mist  be  tried,  and  the  threshold  value  for 
the  nuiher  of  carpatible  point  necessary  to  inply  that  the  solution  is  valid.  Wfe  can  set  the 
first  at  10°  as  determined  fran  the  expected  error  for  a  10  microsecond  average.  The  third 
parameter  should  be  greater  than  10  to  ensure  that  the  equivalent  100  microsecond  average  is 
obtained  The  second  parameter  would  have  to  be  determined  at  prototype  stage  by  the  available 
processing  power  and  actual  prototype  tests  on  data  to  indicate  performance. 

In  suimary,  a  single,  hal  f  wave  I  ength  interferometer  sys  ton  should  provide  angular  resolu¬ 
tion  of  2°  with  a  time  resolution  of  about  one  hal  f  mi  I  I  i  second  For  a  typical  flash  of  one  hal  f 
second  we  wou I d  have  one  thousand  data  points.  From  consideration  of  the  duty  cycle  for  radia¬ 
tion  given  by  hhyenga  and  Warwick  (1981)  and  Rhoded  (1985)  we  expect  some  20  to  100  to  pass  the 
consensus  algoritlmas  valid  data  to  be  passed  on  for  display. 

2.  Configuration  and  Angular  Coverage.  Possible  interferometer  configurations  and  the 
resultant  angular  coverage  are  very  difficult  to  specify  in  any  general  fashion  Allowable 
antenna  type  and  location  points  will  vary  considerably  from  one  aircraft  to  another.  The  site 
errors  and  potential  shadowing  fran  the  craft  will  then  also  be  craft-dependent.  A  horizontal 
baseline  for  azinuthal  coverage  is  easier  to  find  on  a  craft  than  is  a  vertical  baseline  for 
elevation  coverage.  For  the  sake  of  this  feasibility  study  we  present  the  sinplest  configura¬ 
tion,  discuss  the  coverage  and  errors  that  result  and  suggest  additions  which  could  enhance 
performance. 

We  assume  a  source  can  arrive  from  any  azinuth  but  will  be  I  imited  to  about  +  20°  in 
elevation  during  normal  operation  as  shown  in  Appendix  I.  The  sirrplest  configuration  would  be 
two  antennas  separated  by  one-half  wavelength  located  on  the  leading  edge  of  the  wings  of  the 
craft,  thus  limiting  the  patterns  to  the  forward  direction  (Johnson  and  Jasik,  1984,  chapter  37). 
This  configuration  would  not  be  capable  of  elevation  angle  determination. 

Figure  12  shows  lines  of  constant  phase  measurement  for  the  forwad  180°  in  azinuth  and  +  20° 
in  elevation.  The  measured  phase  equals  the  product  of  the  sine  of  the  azinuth  angle  and  the 
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cosine  of  the  elevation  angle  (Hayenga  and  Warwick,  1981).  Sources  not  at  0°  azlnuth  with 
respect  to  the  aircraft  would  lead  to  an  error  in  azimuth  determination.  For  exzrple,  consider  a 
source  at  -60°  azinuth  at  60  ton  range,  8  ton  above  ground  received  by  a  craft  at  15000  ft.  From 
Appendix  I,  the  direct  wave  would  arrive  at  +4°  el  and  the  ground  wave  at  -12°  el.  Figure  shows 
the  sources  and  the  possible  source  azinuth  determinations  which  vwuld  be  made  from  each  wave  and 
the  the  corrtbined  wave.  The  direct  wave  would  be  located  at  57.6°,  the  ground  wave  at  57.9°,  and 
carbined  wave  at  58.7°.  The  errors  due  to  geometry  In  this  exarrple  are  only  slightly  larger  than 
the  expected  system  errors.  At  azinuth  angles  less  than  +  60°  the  errors  decrease.  For  the  sane 
configuration  except  on  azinuth  of  80°  the  sources  are  fixed  at  70°,  74°,  76.5°;  and  the  errors 
are  now  significant. 

It  should  be  noted  frcm  Figures  12  and  13  that  the  largest  geometry  errors  occur  at  large 
azinuth  and  elevation  angles.  In  other  words  errors  increase  outward  fron  the  flight  path.  The 
larger  errors  at  the  edges  of  the  geometry  would  not  be  as  critical  as  if  they  had  occurred  near 
the  f  I  ight  path. 

In  stannary,  a  single  horizontal  base  line  should  be  capable  of  providing  azinuth  angles 
within  the  2°  expected  error  of  the  syston  over  ♦_  60°  in  azinuth  fron  the  flight  path.  This 
anticipates  that  potential  site  errors  can  be  reduced  to  the  same  level  during  prototype  develop¬ 
ment.  Consideration  of  a  crossed  baseline  (i.e.,  vertical)  for  elevation  determination  can  be 
done  from  the  previous  discussion  by  rotating  the  diagrams  by  90°. 

C  Range  Determination. 

In  their  survey,  Parker  and  Kasemir  (1980)  suggest  that  two  interferometers  placed  at  the 
wingtips  of  an  aircraft  might  be  used  to  determine  range  by  triangul atiorv  If  r  is  the  spacing 
of  the  two  systems,  to  first  order  the  angular  resolution  needed  to  triangulate  at  range  R  is 
given  by  the  angle  that  r  subtends  at  range  R.  In  degrees  the  angular  resolution,^,  for 
triangulation  is  approximately 

A*  -  5.73  J  . 
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It  is  easy  to  show  that  triangulation  is  not  feasible.  If  we  assure  a  moderate-sized  craft  with 
r  =  30  meters,  the  2°  angular  resolution  quoted  in  the  previous  azinuth  determination  section 
would  allow  ranging  to  1  knv  To  range  to  ISO  km  would  require  an  angular  resolution  of  j01°. 
This  is  clearly  not  feasible  by  any  method  know)  to  the  author. 

In  the  phenomenology  section  we  presented  results  from  studies  of  radar  returns  frcm  light¬ 
ning.  By  scaling  the  antenna  size  frcm  those  studies  dowt  to  size  suitable  fran  aircraft  and 
than  scaling  up  the  required  transmitted  power  to  accomodate  the  reduction  in  antenna  gain  and 
keeping  all  other  factors  constant,  a  radar  of  1-10  megawatts  pulse  power  would  be  required  to 
see  lightning  channels  out  to  150  tars  We  can  conclude  that  radar  is  not  feasible  for  ranging  on 
lightning. 

Amethod  for  range  determination  using  the  radiation  from  I  ightning  which  might  be  feasi¬ 
ble  uses  the  path  length  difference  in  propagation  of  the  direct  and  reflected  waves  frcm  the 
source.  This  path  difference  appears  at  the  aircraft  as  a  time  difference  between  the  direct  and 
reflected  waves  which  would  be  identical  signals  except  for  changes  upon  reflection  of  the 
reflected  wave.  Studies  presented  by  Reed  and  Russell  (1953)  showed  little  change  in  pulse 
shapes  between  direct  and  reflected  signals.  Markson  (personal  comunication,  1985)  is  studying 
the  effects  of  reflection  on  lightning  signals. 

typendix  II  shows  plots  of  the  path  difference  and  time  delay  between  the  waves  for  the 
system  geanetry.  The  difference  is  a  function  of  source  height,  source  range,  and  aircraft 
altitude.  Since  the  aircraft  altitude  is  known,  range  determination  invo  I  ves  measurement  of  the 
time  delay  and  source  height  or  elevation  angle  seen  at  the  craft. 

Measurement  of  the  time  delay  is  very  straightforward.  The  inpulsive  lightning  signal  above 
would  not  have  a  strong  autocorrelation  over  the  one  to  forty  microseconds  of  the  the  delay  times 
given.  If  we  perform  an  auto-correlation  on  the  received  signal  vtftich  consists  of  the  direct  and 
ground  reflected  signals,  we  can  expect  the  strongest  autocorrelation  peak  to  occur  at  the  path 
delay.  In  effect,  the  direct  and  ground-reflected  signals  are  being  cross-correlated.  The 
inp  I  (mentation  of  the  autocorrelation  could  be  done  with  analog,  digital,  or  hybrid  technology  in 


the  prototype  stage.  As  this  is  a  rapidly  advancing  technology,  it  would  be  premature  to  suggest 
a  technology.  Fran  the  plots  we  conclude  that  a  resolution  of  about  200  ns  in  the  first  10 
microseconds  and  1  microsecond  fran  10  to  40  microseconds  \K>uld  be  sufficient  for  the  correlator. 

This  range  determination  method  critical  ly  depends  on  the  height  of  the  source.  For  exarrple 
consider  a  measured  delay  of  2  microseconds  when  the  aircraft  is  at  15,000  ft.  At  the  extremes, 
this  could  be  due  to  a  source  2  km  above  ground  at  30  km  or  a  source  10  bn  above  ground  at  150  bn 
as  shown  in  Figure  A-ll-3.  Using  terrain  information  and  the  phenaneno logical  consideration  that 
VtF  lightning  sources  occur  mostly  fran  8-12  km  MSL  we  might  reduce  the  height  spread  to  four 
kilometers,  say  4-8  km  above  ground.  The  possible  range  to  the  source  is  now  fran  GO  to  120  tou 
Figure  l-A-3  shows  very  little  change  in  elevation  angle  of  arrival  for  sources  from  4  to  8  bn 
above  ground  in  the  range  from  60  to  120  tars  It  does  not  appear  feasible  to  determine  range  at 
long  distances  by  this  method  since  there  is  no  accurate  estimate  of  the  elevation  angle  of 
arrival . 

At  closer  range,  where  accurate  range  is  more  important,  the  situation  is  nuch  better.  As 
the  plots  in  Appendix  I  dononstrate  the  elevation  angles  of  arrival  for  various  sources  heights 
spread  out  at  shorter  ranges.  This  allows  a  better  estimate  of  source  height  for  a  given 
angular  resolution  of  the  elevation  angle  of  arrival.  As  noted  in  the  previous  section,  a  sourc. 
will  usually  appear  at  three  elevations,  frcm  the  direct,  reflected,  and  contained  waves.  These 
estimates  of  source  height  together  with  the  time  delay  would  be  sufficient  to  dev  lop  an  esti¬ 
mate  of  the  source  range.  For  exanple,  consider  a  measured  time  delay  of  9  microseconds  with  the 
aircraft  at  15,000  ft.  Source  heights  and  range  vary  fran  a  4  bn  high  source  at  15  tan  range  to  a 
10  km  high  soruce  at  30  km  range  (Appendix  II  Figures).  The  10  km  high  source  direct  wave  would 
arrive  at  ♦  11°  elevation  whereas  the  4  bn  high  direct  wave  would  arrive  a  -3°  elevation,  a 
difference  easi  ly  measured  (Appendix  I  Figures). 

The  above  exanples  illustrate  that  the  ground  reflection  technique  is  feasible  for  close 
ranges  (<-60  tan)  but  not  very  accurate  at  longer  ranges.  If  the  technique  is  inplemented  in  a 
prototype  design  a  ccnplete  error  analysis  would  be  necessary. 


It  is  irrp licit  that  this  technique  requires  the  measurement  of  elevation  angles  of  arrival 
using  a  vertically  oriented  interferaneter  baseline.  Vfe  would  recomnend  this  for  a  prototype 
stage.  However,  research  on  the  location  of  VFF  sources  fron  lightning  indicate  that  al  I  flashes 
may  originate  near  the  negative  charge  layer  which  has  also  been  tied  to  the  envirormental 
terrperature.  If  these  indications  can  be  verified  by  further  basic  research,  it  would  be  feasi¬ 
ble  to  imke  a  range  determination  using  only  the  initial  few  milliseconds  of  a  flash  and  an 
assured  height.  This  could  eliminate  the  need  for  the  vertical  baseline  and  its  attendant 
doubling  of  system  ccrrplexity. 

In  suimary,  the  ground  reflection  technique  for  range  determination  fa  I  Is  off  in  accuracy  at 
increasing  range.  The  technique  is  untested  in  that  no  actual  data  of  lightning  radiation 
received  at  an  aircraft  with  reflected  and  direct  waves  exists.  Tests  should  be  performed  on  the 
effect  of  reflection  on  lightning  signals  in  the  lower  VKF  range  and  the  correlation  with  the 
direct  wave  as  a  first  step  in  prototype  developrent. 


V.  Conclusions 


The  feasibility  of  an  airborne  lightning  location  system  using  the  VFF  radiation  from  light¬ 
ning  has  been  examined  in  this  report  in  three  areas.  Fran  the  phenomenological  standpoint  we 
have  shewn  that  our  current  understanding  of  the  WF  radiation  fran  lightning  is  still  very 
limited  but  appears  premising  for  a  detection  system  Mb  have  shawm  that  there  is  little  estab¬ 
lished  connection  between  lightning  location  and  other  potential  weather  hazards  from  thunder¬ 
storm  such  as  turbulence  or  hail.  Further  basic  research  on  thunderstorms  oust  be  performed 
before  we  can  use  lightning  location  as  an  indicator  for  other  hazards.  The  pherxmenical  con¬ 
siderations  of  source  size,  movement  and  variability  in  extent  fran  flash  to 
flash  are  basic  limitations  to  any  syston  No  new  information  can  be  gained 
by  a  more  technical  ly  advanced  system 

There  are  significant  potential  problems  which  nxist  be  considered  as  the  lightning  signal 
propagates  to  the  aircraft.  The  syston  mist  be  considered  as  a  complete  link  including  source 
characteristics,  propagation,  and  receiver  design.  Again,  the  lack  of  basic  research  on  VFF 
lightning  characteristics  limits  the  information  available  for  syston  design  sanewhat,  requiring 
sane  arpirical  work  at  a  prototype  stage. 

Fran  the  technical  standpoint  we  have  shawm  that  an  interferometric  syston  operating  in  the 
30-100  Mfz  range  with  1  NHz  bandwidth  is  feasible  for  locating  lightning  in  the  range  fran  10  to 
150  km  Angular  coverage  of  +  60°  with  2°  resolution  appears  feasible  with  a  single  baseline  of 
one  half  wavelength.  A  second  baseline  would  resolve  elevation  if  desired  Range  determination 
by  ground  reflection  is  feasible  with  resolution  -5  Ion  for  ranges  less  than  ~60  Ion  but  degrades 
severely  at  longer  ranges.  The  technical  considerations  are  tied  to  our  current  phenaneno logical 
understanding.  It  should  be  stressed  that  current  constraints  are  not  technical,  but  phenomeno¬ 
logical. 

Hardware  specific  problems  were  not  addressed  These  include  antenna  type  and  size  and 
possible  siting  errors  on  the  craft.  It  Is  assured  that  these  can  be  specified  sufficiently  wei  I 
in  a  prototype  to  not  degrade  the  given  performance  expectations.  Tolerances  of  carponents  in 


the  interferometer  were  not  addressed.  From  our  experience  with  the  construction  of  a  ground 
station  these  wi  1 1  not  be  a  problem  in  the  airborne  system 


VI.  Recorrmendat  ions 


Two  recorrmendat  ions  onerge  readily  fran  this  study.  First,  the  study  was  limited  by  our 
current  understanding  of  the  VfF  radiation  fran  lightning  Although  significant  advances  have 
been  made,  more  remains  to  be  done,  especially  in  areas  of  relating  the  sources  of  this  radia¬ 
tion  to  the  thunder s tonn  environment.  Not  only  would  such  continued  research  help  further  define 
the  possibilities  of  an  airborne  systen\  such  research  would  also  aid  in  the  testing  of  an 
airborne  system  We  anticipate  difficulty  in  testing  a  system  which  detects  a  still  sonewhat 
unknown  phenomena  in  the  sense  of  sorting  out  technical  problem  from  phenomenological  varia- 
t ions. 

The  second  reconmendat ion  concerns  the  construction  of  the  prototype  unit.  We  do  not 
consider  theoretical  considerations  sufficient  to  determine  fully  the  feasibility  of  an  airborne 
system  We  have  attempted  to  address  these  aspects  which  affect  operation  of  this  system,  but 
can  obviously  nake  no  claim  about  having  seen  every  possible  problem  The  construction  and 
testing  of  a  prototype  system  could  more  accurately  address  a  mrrtoer  of  areas  ehich  were  only 
addressed  as  possible  problsns  or  areas  needing  more  work. 

Specif ical  ly,  we  recoimend  a  prototype  crossed  baseline  interferometer  system  with  one  half 
wavelength  baselines  operating  in  the  region  frcm  30-100  Miz  which  has  tolerable  interference 
levels.  This  wou  Id  a  I  low  assessment  of  azinuth  and  elevation  determinations. 

We  suggest  a  bandwidth  of  1  MHz  with  integration  by  the  methods  discussed  of  half  milli¬ 
second  intervals.  We  also  reconmend  studies  on  the  ground  ref  I  ected  waves  from  lightning  as  seen 
fran  an  airborne  receiver  as  first  test  of  the  ranging  method  described.  As  a  second  step  a 
correlating  receiver  would  be  integrated  with  the  interferaneter  with  appropriate  pro-cessing  to 
determine  range  fran  elevation  and  time  delay  information.  We  also  stress  the  need  to  study 
carefully  the  means  by  vrttich  the  performance  of  such  a  prototype  sys ton  would  be  evaluated. 
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/Appendix  I 


Elevation  Angles  In  Systen  Geanetry 

This  appendix  contains  six  plots  of  the  angle  of  arrival  In  elevation  of  sources  within  the 
system  geanetry.  The  software  used  to  generate  the  plots  is  also  listed  For  the  systsn  geome¬ 
try  we  assumed  the  distance  to  the  lightning  source  would  range  fran  10  ion  to  ISO  km  Source 
heights  were  assured  to  range  fran  2  bn  to  10  Ion  above  ground  with  the  plots  showing  source 
heights  in  two  kilometer  Incronents.  Each  plot  is  for  a  fixed  aircraft  altitude  ranging  from 
5000  ft  above  ground  to  30,000  ft  in  5000  ft  incronents. 

A  flat-earth  geanetry  was  assured  for  this  study  and  the  direct  and  specularly  reflected 
waves  are  shown  In  actual  design  a  spherical  earth  trust  be  used  These  plots  also  have  no 
corrections  for  refraction  or  anomalous  propagation.  Such  factors  rrust  also  be  treated  carv- 
pletely  in  actual  syston  design  The  effects  vmuld  not  significantly  affect  conclusions  derived 


from  these  figures 
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Elevation  angles  of  arrival  for  sources  with  craft  at  15,000  ft. 
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PROGRAM  BOUNCE 

THIS  PROGRAM  CALCULATES  THE  ANGLES  OF  ARRIVAL  OF  THE 
DIRECT  ANO  REFLECTED  WAVES  FROM  AN  ELEVATED  RADIO  SOURCE 
TO  AN  AIRBORNE  RECEIVER. 
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PRINTER  IS  CRT 
PRINT  CHR$< 12  ) 


THE  FOLLOWING  LINES  OF  CODE 
GENERATE  AND  PLOT  THE  AN6LES  OF 
ARRIVAL  FROM  THE  SOURCE  OF  THE 
DIRECT  AND  GROUND-REFLECTED  WAVES. 

in :  1 

GCLEAR 

INPUT  "ENTER  AIRCRAFT  ALTITUDE  IN  FEET" .Alt 
H2=A1 t  * . 3048 
VIEWPORT  10,100,10,90 
Res*-"LOW" 

INPUT  "ENTER  Y  AXES  RESOLUTION  DEF-LOW"  .Res* 
IF  Res$»"LOU"  THEN 

WINDOW  0,1 .S5E+5.-9I  .91 
AXES  15000,10,0.0 
GOSUB  Low 

ELSE 

WINDOW  0,1 .55E+5.-10.5.10.5 
AXES  15000,1  .0,0 
GOSUB  High 
END  IF 
LINE  TYPE  1 
FOR  Y-2  TO  10  STEP  2 
FOR  X-10  TO  150 
D=X  * 1 000 
H  1  =  Y*  1  000 

Theta— (ATN<  (H2-H1  >/D  >  ) 

MOVE  0,0 
PLOT  D. Theta  .1 
LORG  5 
CSIZE  4 

IF  (X  MOD  75  )  =  0  THEN  LABEL  Y 
NEXT  X 
NEXT  y 

LINE  TYPE  1 
FOR  Y-2  TO  10  ?TEP  2 
FOR  X-10  TO  150 
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550 

HI -Y* 1000 

560 

Q*X« 1 000 

570 

Phi=ATN< (H2+H1  I/O ) 

580 

PLOT  0  Phi ) , 1 

590 

L0R6  5 

600 

CSIZE  4 

G  1 0 

IF  (X  MOO  701=0  THEN  LABEL 

620 

NEXT  X 

630 

PENUP 

640 

NEXT  Y 

650  ! 

660 

60SUB  X_axes 

670 

INPUT  "ANY  FURTHER  PLOTS?"  .Answer* 

680 

IF  UPC*< Answer* )=”Y"  THEN  Again 

690 

STOP 

700  Low: 

i 

710  ' 

720  ' 

730 

VIEWPORT  0,120,10,90 

740 

WINDOW  0,120.-91 ,91 

750 

FOR  I*>-9  TO  9 

760 

MOVE  7,I*'0 

770 

L0R6  5 

780 

CSIZE  2.SJ 

790 

LABEL  1*10 

800 

NEXT  I 

810 

MOVE  4,0 

820 

CSIZE  4 

830 

LOIR  90 

840 

LABEL  ’ANGEL  OF  ARRIVAL  < 0E6  >“ 

950 

LDIR  0 

860 

VIEWPORT  10,100,10,5® 

870 

WINDOW  0.1 .55E+5.-91  ,91 

880 

RETURN 

890  High 

:  i 

900  ! 

910  i 

920 

VIEWPORT  0,120,10,90 

930 

WINDOW  0,120,-10.5,10.5 

940 

FOR  I=-I0  TO  10 

950 

MOVE  7.1 

960 

CSIZE  2.5,1 

970 

L0R6  5 

980 

LABEL  I 

990 

NEXT  I 

1000 

CSIZE  4 

1010 

LOIR  90 

1020 

MOVE  2,0 

1030 

LABEL  " AN6LE  OF  ARRIVAL  <OEG>" 

1040 

LOIR  0 

1050 

VIEWPORT  10,100,10,90 

1060 

WINOOW  0,1 .55E+5.-10.5.I0.5 

1070 

RE  TURN 

1080  X_axes:  1 
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VIEWPORT  0,120,0.100 
WINDOW  0,120,0,100 
MOVE  1 1 1  .50 
CSIZE  3.0 
LABEL  “RANGE  (Km)" 
VIEWPORT  10,100,0,10® 
WINDOW  0,1 . 55E+5  ,0 , 1 00 
FOR  X-10  TO  150 
I=X*1000 

IF  (X  MOD  30  )*0  THEN 
MOVE  I  ,48 
LABEL  X 
END  IF 
NEXT  X 

VIEWPORT  0,120,0,100 
WINDOW  0,120,0,100 
MOVE  80,80 
LABEL  “DIRECT  WAVE" 

MOVE  80,10 

LABEL  “REFELCTEO  WAVE’ 
MOVE  70,90 
CSIZE  4.5 

LABEL  “AIRCRAFT  ALTITUDE" 

RETURN 

ENO 


; Alt ;  “Ft 


Appendix  II 


Path  Differences  in  System  Geometry 

This  appendix  contains  twelve  plots  shown  as  six  figures  presenting  the  path  length  differ¬ 
ence  or  time  delay  for  reception  of  a  signal  in  the  system  geometry  propagating  directly  and 
propagating  by  means  of  a  single,  specular  reflection.  The  software  used  to  generate  the  plots 
is  also  listed.  The  systan  geanetry  and  assurptions  are  identical  to  those  given  in  Appendix  I. 
The  lower  plot  in  each  figure  is  a  plot  of  the  shorter  differences  (longer  range  sources)  expand¬ 
ed  for  clarity. 
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Figure  II-A.l.  Path  difference  and  time  delay  for  craft  at  5000  ft. 
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Figure  IT-A.3.  Path  difference  and  time  delay  for  craft  at  15,000  ft, 
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Figure  II-A.U.  Path  difference  and  the  delay  for  craft  at  20,000  ft. 
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PATH  DIFFERENCE  (Km>  PATH  DIFFERENCE  (Km  J 


10 

PROGRAM  PATH_QIr 

20 

This  PROGRAM  CALCULATES  THE  DIFFERENCE 

INPATH  LENGTH 

30  ' 

BETWEEN  A  DIRECT  WAVE  AND  REFLECTED  WAVE,  OR  THE 

40 

PATH  DIFFERENCE  BETWEEN  INCOMING  WAVES 

TO  AN  AIRBORNE  RECEIVER 

50 

r*  -x 

VERSION  1.02.  C. Rhodes.. 

70 

0DT I ON  BASE  > 

w  V 

DIM  Path_di f <  S ,30  • 

90 

PLOTTER  IS  3," INTERNAL" 

100 

GRAPHICS  ON 

1  '0 

GCLEAR 

120 

THE  FOLLOWING  LINES  OF  CODE  ARE  USED 

t  7  ^  1 

J  V 

TD  ENTER  THE  VARIOUS  PARAMETERS  FOR 

*  •  7» 

“*  V 

CALCULATING  THE  PATH  DIFFERENCE 

•  C.  ■* 

1  w  V 

BETWEEN  THE  DIRECT  AND  REFLECTED  WAVE. 

1  50 

•  "Irl 

T~E  c  LAT  EARTH  APPROXIMATION  IS  USED. 

<  X.’ 

i  e0  Aga  j 

r.  •  1 

190 

GCLEAR 

230 

INPUT  ENTER  THE  RECEIVER  ALTITUDE  IN 

FEET" ,Rc  vr_a 1 1 

210 

Rc  vr_a 1 1  *Rc vr_a It* .3048IM/Ft 

223 

*2=Rcvr_alt 

230 

FOR  V=2_TC  10  STEF  2 'KM 

240 

EOF  a =10  TO  1 50  STEP  5 'KM 

253 

I  =  V  /  2 

253 

J  =  X  /  5  -  1 

270 

y | = y* ! 

2  f0 

X 1 =X  *  1 000 

293 

R1 =SQRi i V2-Y1  )  2  +  X  1  *X  1  ) 

~.‘V* 

w’  X.  X^ 

R2  =  SQR< < Y2  +  Y1  )  ' 2  +  X 1  * X 1  ) 

7  :  i“>i 

V 

Fat H_di f  < I ,J  )=R2-R 1 

::o 

NEXT  X 

330 

NE.xT  t 

3-‘*0 

M3a_J1 f -0 

*50 

FOR  1=1  TO  5 

.1  z  C 

FOR  J=!  TO  29 

7”  1 

.  X. 

Ma-_d;f=MAV( ABS<Path_dif( I ,J  >  ) 

,Ma - _di f i 

T  ** 

»'Ct 

r :  E  V  t 

393 

NEXT  I 

-100 

»i  ‘  0  1 

THE  FOLLOWING  LINES  Oc  CODE  SCALE  AND 

430 

Plot  THE  RESULTING  PATH  DIFFERENCES 

1  7  "X 

VERSUS  RANGE  FOR  THE  VARIOUS  SOURCE 

440 

ALTITUDES 

INPUT  "THE  Y_h*es  TICK  RESOLUTION"  .Tick s 

AXES  15000 ,TicK s  ,0  ,0 

Micro_res=<  Y_sc  a  1  e  ■■  3  . E  +  9  '  •'(  !  .  0E-6  > 

Micro_re5 1 -( i _?cal e I  '3 . 0E  +  9  '  ■  1  . 0E - B  > 

WINOOW  0  ,  !  . 5EE  +  5  ,0  ,M l cro_r es I 

AXES  0  ,  .  1  *M i c ro_re 5  ,  I  . SEE * 5  , 0 

WINDOW  0  , ! . 55E+5  ,0  ,  Y_sc a  1 e I 

FOR  I = !  TO  5 

FOP  J  =  1  TO  29 

LINE  TYPE  I+4 

PLOT  i  J  + 1  )*5*  I  000  .Pa  t  h_d  l  f t  I  ,  J  ) 
NEXT  J 
MOVE  0,0 
NEXT  I 

LABEL  THE  X-AxES  IN  Kn 

LORG  5 
LDIR  0 
CSIZE  3.5 
LINE  TYPE  I 
VIEWPORT  1  2  .  H0  ,0 , 100 
WINOOW  0,1  .SSE  +  5, 0,100 
FOF  1=0  TO  10 
X  ■.  =  1  *  1  5000 
X  x =  I  *  1  5 
MOVE  X*  ,5 

LABEL  USING  "ODD"iXxx 
NEXT  I 
BEEP 

VIEWPORT  0,120,0,100 
WINDOW  0 . 1 20  .0  ,  '00 
MOVE  50  .3 

LABEL  "PANGE  FROM  SOURCE  (Kn)" 

LABEL  THE  Y-AXES  IN  Kn 

VIEWPORT  0 , 1 20 , 1 0  ,95 

WINDOW  0 , 1 20  ,0 , Y_sca 1 e 1 

FOP  1  =  1  TO  INK Y_3caie/TicKs  : 

Y/= I • T  i c(  5 

Y/>  =  I  * T i cl  a /  I  000 

MOVE  ?.5,i'/ 

LABEL  Yyy 
NEXT  I 
LDIR  PI/2 
VIEWPORT  0,120  ,0  , 

WINOOW  0 , 1 20  ,0  ,  I  00 
MOVE  2.0,50 

LABEL  "PATH  DIFFERENCE  '.Kr>“ 

LABEL  SECOND  /-AXES  IN  MICROSECOND 


■UWIf^rJ  ■iwvr'l 


vmr 


1  060 

VIEWPORT  0,1 20  .'0,55 

1090 

WINDOW  0 , 1 20  ,0  ,Mi cro_res 

1  100 

LOIR  0 

1  1  10 

FOR  1=1  TO  "3 

l  1  20 

Y,=e 

'  130 

>  >■=  I  «  .  1  *M icro_^s5  ) 

i  ’  4? 

MOVE  112,  Yy 

1  150 

L0R6  2 

1  150 

CSIZE  3.0, .6 

;  ril 

1  KJ 

LABEL  USING  " D . D ” ; ' 

1  1  20 

NE*T  I 

l  I  90 

WINDOW  0,100,0,100 

’  200 

LORG  5 

’"10 

MOVE  119,50 

1  220 

LDIR  ° I / Z 

1230 

CSIZE  3.5,-" 

1  240 

LABEL  “TINE  DIFFERENCE  (i 

’  250  I 

1 250  1 

GENERATE  GRAPH  LEGEND 

1270  i 

i  250 

line  type  > 

i  2  50 

LDIR  0 

l  200 

VIEWPORT  50,90,50,90 

1  3 ' 0 

FRAME 

1  320 

CSIZE  3.0 

'330 

VIEW?GPT  0,120,0,  !00 

,  -40 

WINDOW  0,120,0,130 

'350 

MOVE  70 ,9S 

1  350 

LCPG  5 

i  3 •'.y 

LABEL  “SOURCE  HEIGHT  km" 

1350 

LORG  1 

i  ~S>~* 

_  .  V. 

MOVE  53. S0 

’  400 

CSIZE  4.5 

M>0 

LABEL  “2Kh  = " 

i  a  2  0 

LABEL  “ 4Kr 

'430 

LAEEL  "GKr  = “ 

’  440 

L  "  B  £  L  “  S  K  r  =  " 

,  4C? 

■  «  5  E  L  “  1  0  *■  =  “ 

<  '50 

FOR  1=1  TO  5 

'  4-0 

MOVE  £5,61-'  I-t  '*4 

:  44  0 

LINE  TYPE  1+4 

’4  90 

C0R  2=55  70  99 

1500 

DRAW  2  ,5  '  1-1  )* 

it  ■(> 

NE  *  T  Z 

1  520 

FENUP 

1520 

V£*T  I 

1  540 

1550 

’.'2'W  LABEL  .  hE  hRAPH 

i  550 

i 

i  c-^ 

w  V. 

VIEWPORT  0 , 1 20 ,0  ,  '03 

!  5E? 

WINDOW  0 , 1 20  ,0 . 1 00 

t  c 

MOVE  60,95 

'£00 

LORG  5 

'510 

_ I NE  T  > PE  ' 
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LDIP  0 
CSI2E  4.0 

LABEL  '  AIRCRAFT  ALTITUDE"  “  ;Pcvr_alt ;  "rt 


THE  FOLLOWING  LINES  OF  CODE  DEFINE  THE 
PLOTTING  DEVICE .  IN  GENERAL  THIS  WILL 
BE  THE  THINKJET  PRINTER  AT  HF-IB 
ADDRESS  AT  70S 


DUMP  DEVICE  IS  706  .EXPANDED 
PRINTER  IS  705 

INPUT  "GENERATE  A  HARO  COP  <  PLOT?”  .Answer* 
IF  UPC*< Answer*  :  =  "Y“  THEN 
FOP  I"!  TO  8 
PRINT 
NEXT  I 

DUMP  GRAPHICS 
PRINT  CHR*( 12 ) 

END  IF 

PRINTER  IS  CRT 
BEEP 

INPUT  "ANY  FURTHER  ALTITUDES?"  .Answer* 

IF  UPCSI Answer* >»"Y*  THEN  GOTO  AQain 
STOP 
END 


